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We study theoretically the contribution of fluctuating Cooper pairs to the persistent current 
in superconducting rings threaded by a magnetic flux. For sufficiently small rings, in which the 
coherence length £ exceeds the radius R, mean field theory predicts a full reduction of the transition 
temperature to zero near half-integer flux. We find that nevertheless a very large current is expected 
to persist in the ring as a consequence of Cooper pair fluctuations that do not condense. For 
larger rings with R ^> £ we calculate analytically the susceptibility in the critical region of strong 
fluctuations and show that it reflects competition of two interacting complex order parameters. 

PACS numbers: 74.78.Na, 73.23.Ra, 74.25. Ha 



Introduction and main results- Superconducting fluc- 
tuations have been the subject of intense research during 
the last decades [l[ . At temperature above the transition 
temperature T c to the superconducting state, when the 
system is still metallic, pairs of electrons are formed for 
a limited time. These superconduting fluctuations affect 
both transport and thermodynamic properties. 

In bulk superconductors T c can be reduced or even 
completely suppressed by various phase-breaking mecha- 
nisms, for example by applying a magnetic field or intro- 
ducing magnetic impurities. A special situation occurs 
for superconducting rings and cylinders threaded by a 
magnetic flux <f>. T c is periodically reduced as a function 
of 4>, a phenomenon known as Little-Parks oscillations Q . 
The period of the oscillations is equal to 1 as a function 
of the reduced flux <p — 4>/<j>Q, where the superconducting 
flux quantum is cj>o = 2irhc/2e = ir/e [3J, see Fig. [TJ 

The magnitude of the maximal reduction in T c is size- 
dependent. As we see in Fig. [TJ mean field (MF) theory 
predicts that for small rings or cylinders with r = R/£ < 
0.6 the transition temperature is equal to zero in a finite 
interval close to half-integer flux, giving rise to a flux- 
tuned quantum phase transition, see also Eq. (jSJl below. 
In this Letter we show that the pair fluctuations give a 
large contribution to the persistent current (PC) I even 
at fluxes for which T c is reduced to zero and the system 
has a finite resistance. 

Recent experiments added significantly to our under- 
standing of fluctuation phenomena in superconductors 
with doubly-connected geometry. Strong Little-Parks os- 
cillations in the region where £ > R, where T c is re- 
duced to zero, have been observed in a transport mea- 
surement on superconducting cylinders [ij]. Koshnick et 
al. Q measured the PC in small superconducting rings 
in the regime where R > £, for the smallest rings under 
study T c was reduced by » 6%. In this Letter we discuss 
the PC both in the regime of moderate T c suppression for 
r = i?/£>las well as the strong Little-Parks oscillations 
for r < 0.6. Before presenting details of our approach, 
we summarize the main results of our analysis. 
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FIG. 1: Mean field (MF) phase diagram. T cv separates the 
metallic (high T) and the superconducting (low T) phase 
as a function of the flux ip = (j>/4>o through the ring. For 
sufficiently small rings with effective radius r = R/£ < 0.6 
MF theory predicts a full reduction of T c for fluxes between 
ifc ~ 0.83r and 1 — <p c near ip — 1/2. The transition line 
reflects the condition — 0, cf Eq. JSJl. 



/. Regime with r = R/£ < 0.6: For r < 0.6 the mean 
field T c vanishes and one would naively expect a small 
normal state PC. We find, however, that close to the crit- 
ical mean field line (see Fig. [JJ there is a parametrically 
large enhancement of the PC due to quantum fluctua- 
tions that decays only slowly away from that line. The 
magnitude for the normal PC is In ~ > where 

g is the dimensionless ring conductance [fj, |7j . Our cal- 
culations show that the PC due to pair fluctuations near 
the critical flux ip c is parametrically larger and at low T 
given by 



JFL « pfog -r— 



(1) 



where Aip = (cp — ip c )/(p c measures the distance to the 
critical flux (p c . When increasing T the PC initially grows 
before going through a maximum at finite T, where it 
can considerably exceed the result of Eq. ([I]) [see Fig. [3]. 
Since r _1 = £/R is a number of order 1 and -j^- = ^^f 
for a weakly disordered superconductor, we find an en- 
hancement factor of log(o) log(l/Ay>). 

Our results are obtained for the case when the flux 
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acts as a pair breaking mechanism. Other pair break- 
ing mechanisms, e.g. magnetic impurities or a magnetic 
field penetrating the ring itself will lead to similar re- 
sults. They cause a reduction of T c to zero, the pair 
fluctuations, however, lead to a parametric enhancement 
of the PC in the normal state. Ref. [1] suggests that a 
similar mechanism due to magnetic impurities is related 
to the unexpectedly large PC in noble metal rings [gL flol| . 

A metallic state with small but finite resistance was 
observed experimentally in superconducting cylinders [H, 
ll| with ip ps 1/2. Further theoretical and experimental 
studies will be needed in order to clarify the relation to 
our findings, where a large PC is caused by pair fluctua- 
tions that are unable to condense. 

77. Regime with r > 1: The case r > 1 is suitable 
for the description of the experiments on persistent cur- 
rents by Koshnick et al. [5(. Previously the experiment 
has been interpreted using a one-dimensional Ginzburg- 
Landau theory to describe the order parameter fluctu- 
ations [13] • Following these lines one has to resort to 
numerical methods 13] in order to describe the critical 
region close to T c , where fluctuations proliferate. 

Our key observation is that part of the rings in the 
experiment allow for a description using a suitable gen- 
eralization [3] of the Od Ginzburg-Landau theory. In- 
deed, following an expansion of the order parameter field 
-0(t9) in terms of angular momentum modes ip n , a simple 
physical picture arises in the limit yfg 3> r. Two of the 
modes compete with each other close to half-integer flux, 
while at the same time both of them strongly fluctuate 
in the critical regime close to T c . 

Formally, the competition arises due to the quartic 
term in the GL functional that induces an interaction 
between the modes [lj| and reveals itself in the experi- 
ment mostly in the "slope" of the PC, the susceptibility 
X = — S • With this insight \ can be calculated analyti- 
cally even in the critical fluctuation regime. 

As an example, denoting the susceptibility at T c and 
zero flux by xo and at T c ^ =1/2 by xi/2, we find 



X1/2/X0 ~ -2.7y/g/r. 



(2) 



Experimentally, a strong enhancement of the magnetic 
susceptibility near p = 1/2 compared to ip » was ob- 
served and Eq. © demonstrates that it is controlled by 
the parameter y/g/r. If it is large, the current will rapidly 
change sign as a function of the flux at half-integer flux, 
leading to a saw-tooth like shape of i v . The full T de- 
pendence of y,„ _ 1 /? is given in Eq. |(7j). For the smallest 
rings in Ref. la UM, yfg w 33r. 

Classical GL functional - After presenting the main 
results in Eqs. (JTJ) and (j^J we now give more details of 
our approach starting with the description of rings with 
only a moderate suppression of T c (i.e. r > 1). 

When the superconducting coherence length £(T) and 
the magnetic penetration depth A(T) are much larger 



than the ring thickness, the system is well described by 
a one-dimensional order parameter field tp (l7j . The 
partition function can be written as a weighted av- 
erage over configurations of the order parameter ip, 
Z = J Dipexp[— T jT\. Introducing angular momentum 
modes as = J2 n ^« e ln ® , where V is the volume 
of the ring, the free energy functional takes the form 
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5 n +k,l+ m ^n^* m ^k^l ■ (3) 



Here we wrote a nip = aT c e ntp , where e ntp = [T — T ntp /T c 
is the reduced temperature and T nip = T c [l — (n~ip) 2 /r 2 ] 
is determined by the sign change of the coefficient a n (ip) 
and can thus loosely be interpreted as the transition tem- 
perature of mode tp n (l8j . The mean field transition oc- 
curs at T cip that is equal to the maximal T n for given 
ip, i.e. at the point where the first mode becomes su- 
perconducting when lowering the temperature (cf. Fig. 
[I}. The Od Ginzburg parameter Gi = {2b/a 2 T c V) 1 / 2 is 
an estimate for the width of the critical regime in the 
variable e n . The parameter yfg/r w \jhr 2 Gi has been 
used when stating our results. Its relevance is now easily 
understood. 1/r 2 is a measure for the typical spacing be- 
tween the transition temperatures T n for different modes, 
since (T - Ti)/T c = (1 - 2p)/r 2 . This spacing should 
be compared to the typical width of the non-Gaussian 
fluctuation region, Gi. If it is large, a theory including 
only one or two angular momentum modes is applicable. 

Persistent current - The persistent current I is found 
from the free energy F = —T\nZ by differentiation I = 
—dF/dcj). The normalized current is given by 

00 _ 

i = I/(TM= Yl -("-^(NAnl 2 ). (4) 



The averaging is performed with respect to the functional 
J- in Eq. ([3]). i„ is periodic in the flux ip with period one. 
Since it is also an odd function of the flux, it vanishes 
when the flux takes integer or half-integer values. 

Case p> w : The most important contribution in 
the regime of non-Gaussian fluctuations close to inte- 
ger fluxes comes from the angular momentum mode 
ip n with the highest transition temperature T nip . One 
may then approximate Eq. J3f by a single-mode and 
calculate with T n = a n \ip n \ 2 + 7^7 | 4 [19j |. This is 
the Od limit of the GL functional [20( where the func- 
tional integral becomes a conventional integral. Indeed, 
performing the integral in polar coordinates, one finds 

Z = ^ ex P( x ™) erfc(x„), where x n = s n /Gi [Hj]. Us- 
ing now Eq. (|H) with one mode only we find 



i n = 4A(n - <p)f(x n ) for p w n. (5) 
Here A = \jr 2 Gi sa 5,/g/r and f(x) - exp ^ g2) - - 



i/7rcrfc(a:) 



21| . We note in passing the high degree of universality 
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implied by this result: All PC measurements will fall on 
the same curve, if the PC - measured in suitable units 
i = 1 1 {T c /<p Q ) - an d the reduced temperature e v = (T — 
T ctp )/T c are scaled as i ^ i-^=, e v — > e v r^fg. The scaling 
function / was given above. This relation is a valuable 
guide in characterizing different rings in experiments. 

Far above T c one obtains as a limiting case the Gaus- 
sian result for a single mode i n « 2(n—tp) /r 2 e nLp , that can 
also be obtained directly by neglecting the quartic term 
in the GL functional. It is known, however, that as soon 
as temperatures are too high, e n ^> 1/r 2 
to sum the contribution of all modes 
one recovers the mean field result imf 



x 



it is important 
22|. Far below T c 
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for the PC in the superconducting regime. An alterna- 
tive route to finding the mean field result would be to 
minimize the full single mode functional, which leads to 



the condition \ip„ 



y/b and then to use Eq. (@}. 



The PC i n in Eq. ([5]) interpolates smoothly between the 
Gaussian and the mean field result. 

Case tp w 1/2; A very interesting situation occurs at 
half integer values of tp. The transition temperatures for 
two modes become equal, their coupling becomes crucial 
(<p ~ 1/2 for definiteness), and we approximate 14 1 

b 



8=0,1 
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4|^o| 2 |^i| 2 



.(6) 



Calculation of the PC in the presence of the coupling re- 
quires a generalization of the approach used for the single 
mode case 14, 23|. In Fig. [5] we display the PC i 2 as cal- 
culated from Eq. ^ for three different temperatures, 
T c i < T < T c , T = T c i and T < T c i. We compare it 
to the MF result as well as to i 2 o obtained by neglecting 
the coupling |i/'o| 2 |'0i| 2 m Eq. {6]). 

Above T cip (e — —0.05 in Fig. [2|), in the region where 
the mean field result vanishes near half-integer flux, the 
PC is purely fluctuational. We deduce from Fig. @ that 
the coupling of the modes is crucial for x(l/2), but not 
for the overall shape when T > T c i . However, just below 
T c i (e = —0.11 in Fig. [2|) the coupling is essential. The 
mean field result is not applicable as it gives an infinitely 
sharp jump in the PC at half-integer flux. The result 
without coupling of the modes, i 2 o, gives a finite slope, 
but it is far from the full current i 2 that includes the mode 
coupling. The coupling drives the current i 2 towards the 
mean field approximation iuF which includes only one 
mode. This occurs because for a repulsive coupling the 
dominant mode suppresses the subdominant one. Indeed, 
if mode n = is dominant then the coupling adds a mass 
term 2&/V r (|')/' 2 |) IV'il 2 to mode n = 1 and reduces its T c . 

Susceptibility - We will now discuss in more detail 
the slope at half-integer flux, which is most sensitive to 
the coupling between the modes below, and to the non- 
Gaussian fluctuations close to T c . Differentiating the ex- 
pression [23I ] for i 2 we obtain 



Gi 

tions <7i (x) 
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24]. The dimensionless smooth func- 
2^e3 a;2 erfc(a;) - 2£ and g 2 (x) = 

2^J(x) - " - 27M e ^ 2 erfc ^) " 1; Wher<3 J ^ = 

f™dt e^' 2 erfc(t), obey 5l (0) w 0.78 and g 2 (0) w 0.315. 
For large A = l/r 2 Gi w h^/g/r one can neglect the first 
term in Eq. Q. Then one obtains X1/2 — ~ 4A 2 g 2 (x). 
For the susceptibility close to integer flux one easily ob- 
tains Xo = 4A/(xo) from Eq. ([5]). Comparing to the 
expression for X\l%i we nn< i Eq. 

This is the strong enhancement of X1/2 compared to 
Xo observed in the experiment J5J . Comparison with the 
numerical calculation of Ref. [5( shows that our analyt- 
ical results are accurate to within a few percent already 
for VsA->8 [HQ. 




FIG. 2: The PC i = I/(T c /<f> ) as a function of the flux ip. 
Parameters are r — R/£ = 1.66, A = l/(r 2 Gi) ~ = 
50, e = (T — T c )/T c . The transition for tp = 1/2 occurs at 
e = —0.09. Full lines: 12 calculated with T of Eq. (|6j), it 
takes into account two modes and the interaction between 
them. We compare 12 to two approximations, which neglect 
this interaction. Dotted lines: The mean field approximation 
imf [discussed before Eq. ©] and dashed line: 120 calculated 
with T of Eq. without coupling [25[. Inset: MF phase 
diagram, superconducting region in grey. 

Quantum critical regime - So far we have discussed the 
limit r = R/£ > 1, where the suppression of T c is small 
and a finite temperature phase transition occurs. We will 
now discuss the case where r — R/£ < 0.6 and T c is re- 
duced to zero at a critical flux ip c near tp = 1/2, see FigO] 
Near the quantum critical point (QCP) it is no longer 
legitimate to use the classical GL functional, in which 
only the static component of the order parameter field 
is considered. Instead, all Matsubara frequencies should 
be taken into account in the imaginary time formalism. 
The full fluctuation propagator is given by: 



("C)nfc=k 



|fi fe |/2 



2ttT 



(8) 



27rfcT are 
Following the 

X v =i/2 — x/(Tc/<f>o) = 4A g\ (x) — 4A^ g 2 (x) , (7) standard approach, we first find the critical line <p(T) 



where a n (ip) = ^T^i( n ~~ v) 2 an d _J?fc = 
bosonic Matsubara frequencies [2" 



27] 
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in the temperature-flux plane by equating C 



oo = 0. 

For the QCP at T — one obtains the critical flux 
ip c = 7rr/(2 v /57i), -f E ~ 1-78 3}. Due to the flux- 
periodicity of the phase diagram, the QCP can only be 
observed in the ring geometry if ip c < 1/2 which implies 
r < ^/27b/7t « 0.6. Notice that this critical value of 
r = R/£ is 20% larger than a naive application of the 
quadratic approximation valid for r> 1 would suggest. 

Restricting ourselves to the interval <p e (0,0.5) we 
find that near the QCP it is sufficient to consider the 
n = mode. In the Gaussian regime we obtain (cf. 
Fig [TJ) the following fluctuation contribution to the PC 
ic = - Tctptfr) T J2k^-ok HI]. Expanding C^ 1 in small 
Aa = [a>o((p) — a((p(T))]/ao(<p c ) we find with logarith- 
mic accuracy i G = -^h(Aa,t) ^ -^j- h (Aa,t), 
h(Aa,t)=]n±; + ±-iP{l + S ),s=^ a iidt = T/T c . 

A few remarks are in order concerning this result. 
The second term in the expression for h is the classi- 
cal Q = contribution to the sum. The upper cut- 
off for the fre quen cy summation has been chosen as 
fl = 2oto(ip(T)) [29(. The function h has the asymptotic 
form h sa j E t/Aa + \n(l/2^ E t) for Aa < t < 1 and 
h « In 1/Aa for t <C Aa < 1. It is important that Aa is 
T-dependent and in order to reveal the full T-dependence 
of %g one should first find the transition line ao (<p(T)). 
h(T) is displayed in Fig. [3] The maximum of h at finite 
T is a result of two competing mechanisms. As T grows 
from zero, thermal fluctuations become stronger. At the 
same time the distance to the critical line becomes larger 
for fixed tp, which eventually leads to a decrease of ic- 
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FIG. 3: The function h that determines the PC close to the 
QCP, iG ~ — ^^r' 1 ! as a function of t = T/T c for given 
Aao ~ 2(93 — tfc) /<fic- <Pc = o — 75= i s th e critical flux at 
T — and 7_e ~ 1.78. h is defined in the text. The dotted 
lines describe classical fluctuations. Inset: MF phase diagram, 
superconducting region in grey. 

Conclusion - In conclusion, we showed that on the nor- 
mal side of the flux-tuned superconductor normal-metal 
transition in small rings the fluctuation PC can be very 
large compared to the normal case and decays only loga- 
rithmically away from the critical point. For larger rings 



as studied in recent experiments we obtained detailed 
analytical predictions for the strong fluctuation region. 
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